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Helium stars as supernova progenitors 



Roni Waldman ^, Lev R. Yungelson ^ and Zalman Barkat ^ 

Abstract. We follow the evolution of helium stars of initial mass (2.2 — 
2.5) M0, and show that they undergo off-center carbon burning, which leaves 
behind '^0.01 Mq of unburnt carbon in the inner part of the core. When the 
carbon-oxygen core grows to Chandrasekhar mass, the amount of left-over car- 
bon is sufficient to ignite thermonuclear runaway. At the moment of explosion, 
the star will possess an envelope of several 0.1 Mq, consisting of He, C, and 
possibl y some H, perhaps pro ducing a kind of peculiar SN. Based on the results 
of Wal dman &: Barkat ()2007D for accreting white dwarfs, we expect to get ther- 
monuclear runaway at a broad range of pc ~ (1 — 6) x lO^gcm"^, depending 
on the amount of residual carbon. We verified the feasibility of this scenario by 
showing that in a close binary system with initial masses (8.5 -I- 7.7) M© and 
initial period of 150 day the primary produces a helium remnant of 2.3 Mq that 
evolves further like the model we considered. 



1. Introduction 

Type la supernovae (SN la) have a relatively small dispersion of luminosity (the 



stand ard deviation in peak blue luminosity is as ~ 0.4—0.5 mag. jBranch &: Miller 



(|1993I )) and are being used as distance indicators ("standard candles"), having 



especial significance in the effort of determining the cosmological parameters of 
our universe. 

The long-standing explanation of the SN la phenomenon has been the ex- 
plosive burning of degenerate carbon in the core of a carbon-oxygen white dwarf, 
which becomes unstable as it grows to Chandrasekhar mass {Mch) either by ac- 
cretion from a binary companion or by a merger of two white dwarfs, following 
the angular momentum loss from the system by gravitational wave radiation. 
However, theoretical models are still far from self-consistently producing an evo- 
lutionary path towards the progenitor and reproducing crucial features of the 
observational data, such a s the compo s ition o f the ejecta. For a deta i led re - 
view of the above see, e.g., iLeibundgutI (|2000l ) : iHillebrandt fc Niemeverl (|2000l ): 



Filippenkd (|2005l ). 



As well, SN la can not be regarded as perfectly homogeneous class, since 
their Hubble diagram exhibits scatter larger than the photometric errors, while 
spectroscopic and photometric peculi arities ha ve been noted with increasing fre- 
quency in well-observed SN la (e.g., iFilippenko 2005). 

Therefore, there is an obvious need for progenitor scenarios that could ex- 
plain the diversity among SN la. Several explanations have been suggested, such 
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as variations in the metallicity of the progenitor, in the carbon to oxygen ratio at 



its center, or in the central density at the tim e of ignition (e.g., iTimmes et al 



20031 : iRopke et al.ll2006l : iLesaffre et aD l2006l ) . The variation of the latter two 



parameters is expected to result from the variation in the initial white dwarf 
mass and in the accretion history. 

In this work we follow the evolution of helium stars with initial mass ~ 
(2.2 — 2.5) Mq and show that they might reach thermonuclear explosion and 
perhaps account for some of the peculiar SNe. 



2. Results 

We followed in detail the evolution of a 2.4 Mq helium star, starting from a 
homogeneous ob j ect. We used the TYCHO evolutionary code described in 
Young Sz ArnettI ( 20051 ). The helium burning convective core has an almost 



constant mass of about 0.82 Mq, and produces a CO core with Xq ~ 0.27. 
Subsequently, a radiative helium burning shell develops above the core, and as 
the core grows to about 1.1 M© carbon is ignited in the core at about 0.3 M© 
(Fig. 1, left panel). After the end of core helium burning the center of the star 
becomes increasingly degenerate, and neutrino cooling is increasingly compet- 
ing with contraction-induced heating. As a result, the maximum temperature 
is attained at an off-center point and it keeps growing until carbon ignites there 
(Fig. 1, right panel). Following off-center C-ignition, the center expands and 
cools, and a series of carbon burning flashes occurs. Eventually burning ceases 
when carbon is almost exhausted in the core, however as can be seen in Fig. 2, 
the innermost 0.5 Mq of the core has a residual Xq ~ 0.02. 

After carbon burning ceases, the center continues to heat and contract while 
the helium burning shell gradually increases the mass of the core. Similarly to 
AGB evolution, the luminosity of the star grows, causing the envelope to ex- 
pand and to develop a deepening convective region, which at certain moment 
penetrates into the helium burning shell and creates conditions for hot bottom 
burning. During this stage the carbon that accumulates below the helium burn- 
ing shell quasi-periodically ignites (see Fig. 1, left panel), similarly to the helium 
flashes occurring in AGB stars. Since very little is known about t he mass loss 
rate of stars of the kind we consider, we applied iReimersI ( 1975l )-based rate, 



which reduced the mass of the star to ~ 1.8 Mq (Fig. 1, left panel). 

Finally, the core grows to Mq^ due to the helium burning shell and carbon 
in the center ignites (Fig. 1, right panel). Convection is initiated, supplying 
more carbon to the central region, and if the amount of carbon is sufficient to 
raise the temperature above the oxygen ignition threshold (~ 1.5 x 10^ K), which 
indeed happens in this case, oxygen will ignite under degenerate conditions and 
initiate an explosion similar to the classical SN la case. 

Models with initial masses (2.2 — 2.5) Mq evolve very similarly to the 2.4 Mq 
case, developing a carbon residue which might ignite thermonuclear runaway. We 
did not follow the evolution further through the hot-bottom burning stage. The 
2.1 Mq model only ignites carbon at a later stage, below the helium burning 
shell, similar to the carbon flashes encountered in the 2.4 Mq case, however, 
since in this case carbon has not been previously depleted in the core, burning 
will probably continue until it reaches the center. We have not followed this 




Figure 1. Evolution of the initially 2.4 Mq helium star model. Left: Kip- 
penhahn plot after the end of core helium burning. Light (red) filled areas 
are convective, dark (blue) line shows maximum thermonuclear energy pro- 
duction rate. Right: Evolution of stellar center on temperature vs. density 
plane (long dash (green) line); also shown is the p — T relation of the inner 
part of the model prior (short dashed (dark blue) and dotted (purple) lines) 
and immediately after off-center carbon ignition (dot-dashed (light blue) line). 
Carbon ignition line for Xq = 0.25 is shown as solid (red) line. 

composition at explosion 




Figure 2. Composition of the initially 2.4 Mq helium star model before 
thermonuclear runaway. 



model further. The models of 2.6 Mq and more massive ones ignite carbon very 
close to or at the center, so that the carbon residue is either non-existent or 
insufficient for thermonuclear runaway. 
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In an earlier work ( Waldman &: BarkatI 20071 ) we explored the similar case 
of CO white dwarfs undergoing off-center carbon burning followed by mass ac- 
cretion. As a function of the amount of residual carbon we got thermonuclear 
runaway at a broad range of pc ~ (1 ~ 6) x lO^g cm~'^. Since the structure of the 
CO core in our helium star models is very similar to that of the mass accreting 
CO white dwarfs, we expect to get a similar result at runaway. 

Supernova which results from thermonuclear runaway in the remnant of 
He-star most probably will not differ photometrically from a "normal" SN la, 
but one may expect presence of strong He-lines in the spectrum, thanks to thick 
He-mantle of pre-SN, making this SN la "peculiar" (N. Chugai, priv. comm.). 

To complete the picture, we tested whether a hydrogen-deficient star very 
similar to our initial models could be created as a result of close binary evolution. 
We begun with a binary of (8.5 + 7.7) Mq, with an orbital period of 150 day. 
After the stage of Roche lobe overfiow, the primary has a total mass of 2.3 Mq, 
a CO core of 1.2 Mq, and a surface hydrogen mass fraction of 0.14. Later, wind 
mass loss and sporadic RLOF reduce the mass to 2.1 Mq. Further evolution of 
the remnant is similar to the above described. 



To summarize, we showed that evolution of helium stars with initial mass about 
(2.2 — 2.5) Mq, which might be produced in close binaries, may suggest a SN 
scenario in which thermonuclear runaway in Mc/^-mass cores is initiated by a 
very small amount of residual carbon, while the stars still have a several 0.1 Mq 
envelope consisting mostly of helium, carbon and possibly some hydrogen. In 
order to give a well justified statement on the observational outcome, our pre- 
SN models should be used for detailed simulations of the explosive runaway and 
spectra modeling. 
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